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ABSTRACT 

We present the first lensing mass measurements of Sunyaev-Zel'dovich (SZ) selected clusters. 
Using optical imaging from the Southern Cosmology Survey (SCS), we present weak lensing 
masses for three clusters selected by their SZ emission in the South Pole Telescope survey 
(SPT). We confirm that the SZ selection procedure is successful in detecting mass concentra- 
tions. We also study the weak lensing signals from 38 optically-selected clusters in ~ 8 square 
degrees of the SCS survey. We fit Navarro, Frenk and White (NFW) profiles and find that the 
SZ clusters have amongst the largest masses, as high as 5 x 10 14 M©. Using the best fit masses 
for all the clusters, we analytically calculate the expected SZ integrated Y parameter, which 
we find to be consistent with the SPT observations. 
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1 INTRODUCTION 

Galaxy clusters are the largest gravitationally bound objects in the 
Universe and can be used as cosmological probes because their for- 
mation and ev olution rate are sensitive to different cosmological pa - 
rameters (e. dEvrardll989i : lHaiman et alJl200ll ; lAIlen et al.ll2004l) . 
The abundance of galaxy clusters as a function of mass N(m, z) 
at high redshift z is particularly sensitive to different cosmologi- 
cal models. To probe cosmology and dark energy we must observe 
galaxy clusters at high redshift and obtain mass estimates for them. 

Observations of the Sunyaev-Zel'dovich (SZ) effect 
( Sunyaev & Zel'dovichl Il98lh are a powerf ul way to prob e 
galaxy clusters by detecting the hot cluster gas (Birkinshaw 19981) . 
SZ-detected clusters are in principle particularly powerful as 
they can be seen to high redshifts. The intensity of the SZ effect 
su mmed over the en tire cluster closely tracks the mass of the clus- 
ter dMotl et alj|2 005). X-ray or SZ effect mass estimates are based 
on simplified assumptions such as a hydrostatic equilibrium for 
the cluster gas. It is becoming increasingly apparent, however, that 
we cannot fully model the complex gas physics in clusters within 
a simple framework. Nonetheless, it will be very challenging to 
calibrate cluster masses at high redshift, so in order to use SZ 
observations to probe cluster properties and cosmological models 
it is important to understand the relationship between mass and 
SZ observables in lower redshift systems. Gravitational lensing 
facilitates the calibration of t he SZ observab l es to o btain accurate 
masses for SZ detections dLewis &Kind 120061: | Sealfo n et al.l 
l200fj) . A large area in the southern sky is currently being surveyed 
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in SZ by the Atacama Cosmology Telescope (ACT) and the South 
Pole Telescope (SPT). 

Gravitational lensing is dependent only on the projected mass 
distribution of the lens and so it is possible to study the mass dis- 
tribution independent of its form, including the distribution of dark 
matter. Gravitational lensing causes small (~ a few %) changes 
to the shape of individual galaxies, which can be used to recon - 
struct the mass distribution in the re gion ^Kaiser & Squires 19931) : 
for a review see lMunshi et al.l d2008h . In this paper we present mea- 
surements of weak lensing masses for clusters which were, for the 
first t ime, detected blind by their SZ decrement JStaniszewski et all 
2008). We also include mass measurements from 24 optically- 
detected clusters. This paper is structured as follows: in Sj2] we 
present the data and discuss the image processing; fj3] describes 
the gravitational lensing methods used; the mass measurements are 
presented in Sj4] and we compare with measurements from other 
techniques and calculate the Y parameters. Throughout this paper 
we assume a flat cosmology with Q m — 0.27, Q,a = 0.73 and 
H = lOOh km s _1 Mpc _1 with h = 0.71. 



2 OBSERVATIONS 

We use publicly-available data from the Blanco Cosmology Survey 
- a National Optical Astronomy Observatory Large Survey Project 
observing 60 nights over 4 years on the Blanco 4m telescope at 
the Cerro Tololo InterAmerican Observatory in Chile. The Mo- 
saic II camera is being used for a deep, four-band optical igriz) 
survey of two 50 deg 2 patches of the southern sky. Two areas of 
southern sky have been targeted, centred on 23 ,l 00 m , -55°12 m 
and 05 h 30 m , -52°47 m . These fields lie within a larger area of 
the southern sky which ACT and SPT plan to survey. The paper is 
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based on observations taken in 2005, with the exception of 2 clus- 
ters from 2006 data. The seeing varies between 0.81" and 1.09" 
with a mean of 0.89". The image reduction was carried out us- 
ing the Rutgers Southern Cosmology Pipeline (flat field correction, 
CCD calibration, removal of saturated star bleed-trails, and bad 
pixel masks). Next the images were aligned, stacked and median 
combined using SWarp (Bertin l2006l ); an astrometric solution was 
found by matching stars to sources in the US Naval Observatory 
Catalog. Additional masks were made to remove saturated stars, 
satellite trails and other blem ishes in the image, remo ving 8% in to- 
tal. For more information see Menante au et alj d2009l) . Note that we 
use AB magnitudes throughout. To calculate photometric re dshifts 
(photo-zs), multi-band SExtractor dBertin & Arnouts| [l996) g,r;i,z 
isophotal magnitudes were used to find redshift probabilit y distri- 
butio ns of each object. This was done u sing the BPZ code dBemtezI 
2000), see also Mena nteau et~ai] d2009l) . We focused on 4 clusters 
found in SZ dStaniszewski et al.ll2008h and also 38 optical clusters 
selected from ~ 8 square degrees JMenanteau et alj|2009l) . Of the 
38 clusters, we found non-zero mass estimates for 24 clusters, and 
upper limits for the remaining clusters. We do not consider 5 of the 
clusters as they are in regions only observed in a single exposure. 



3 WEAK LENSING ANALYSIS 

We used the /-band data, 3 co-added exposures of 450s each, 
for our shear analysis, and measu red galaxy shapes using the 
iKaiser. Squires & BroadhursJ dl995|) (KSB ) method. Our pipeline 
is based on[B acon, Refreg ier & Ellisl |2000) and labelled "MB" in 
iHevmans et al.1 d2006l) . but has been automated to process rapidly 
the large SCS data set. The method deviates from "MB" in auto- 
mated star/galaxy separation and in Point Spread Function (PSF) 
interpolation. We tested the pipeline against s imulated images from 
the Shear TE sting Programme (STEP) (see IHevmans et alj|2006l : 
iMassev et al.ll2007h . Our method underestimates shear, but consis- 
tently throughout a wide range of observing conditions, which is 
henceforth compensated for in our shear measurements by apply- 
ing a calibration factor of 1/(0.82 ± 0.05), similar to "MB". 

As in "MB" our pipeline locates galaxies via the IMCAT0 
HFINDPEAKS algorithm and measures their quadrupole shape mo- 
ments using a Gaussian weight function of width 4r s , where r 9 
is the size of the best-fitting Gaussian. The pipeline fits the shear 
polarizability factor |Tr(P 7 ) as a function of galaxy size. We ex- 
cluded galaxies smaller than 1.1 times the measured seeing, and 
those with S/N < 5, le aving 9 per squar e arcminute to !~23. The 
median redshift is 0.65 dCoil et alj |2004): 69% of galaxies in our 
catalo gue also had photo-zs , which have a consistent z distribu- 
tion dMenanteau et alJ l 2009 V Departing from the "MB" pipeline, 
we perfo rmed star/galaxy s eparation via the automated THELI al- 
gorithm dErben et al. 2005), and automated the removal of galax- 
ies with abnormally large values of shear polarizability, smear po- 
larizability or ellipticity, which was present in the "MB" method 
but labor intensive and slow. The observed shapes of galaxies 
were finally corrected for the blurring effect of the PSF. We mea- 
sured the PSF using the 0.5 unsaturated stars per square arcminute 
with 22.0 > i > 18.1 (S/N 55 to 1670). The PSF ellipticity is 
0.035 ±0.019, where the error is the standard deviation throughout 
the survey. Optical effects and temporal variation of the atmosphere 
and telescope between dithered exposures produce patterns in PSF 
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size and ellipticity. This variation was fitted as a sum of polyno- 
mials (of order 4 in the x and y directions) plus sums of sines and 
cosines (of orders up to 4 in a: and 8 in y). These choices give a 
small r.m.s. residual |ej of 0.0092. 

We estimate masses by fitting an NFW profile 
dNavarro. Frenk & White! 1 19961) to the average tangential shear 
from IMC AT 'S ETPROFILE routine. The concentration index of 
the halo was take n as a function of mass and cluster redshift 
dDolag et alj l2004h . For a robust treatment of missing regions 
of data we use a Wiener-filtered mass reconstruction method 
dHu&Keetonll2002h . To address the mass-sheet degeneracy the 
average surface mass density is set to zero over the whole field 
of view (half a square degree). For the SPT clusters the peaks 
of the surface mass density in the Wiener-filtered maps were, 
where present, assumed to be at the centre of the cluster, otherwise 
(SPT 0547-5345) the luminosity-weighted centre was used. The 
latter procedure was followed for the optically selected clusters. 
Our Wiener-filtered reconstruction was done as follows: the two- 
point correlation function of the lensing convergence, used as prior 
for the Wien er reconstruction, was est i mated from the shear-shear 
correlation dBartelmann & Schneidej 1200 ll) £+(6) in the data 
itself. In order to have a smooth prior, we fitted the measured 
with £+(0) = <x/(l + bff) where a and b are constants determined 
by the fit. To obtain S/N maps of the lensing maps we divided the 
maps by the r.m.s. of 500 noise realisations, which were generated 
by randomly rotating the ellipticities of the sources followed by a 
Wiener reconstruction. The S/N maps were then used to confirm 
the cluster centres found in the mass maps. 



4 RESULTS 

Wiener-filtered mass reconstructions for three of the SPT SZ- 
detected clusters are shown in Figure Q] The fourth SZ clus- 
ter at redshift z ~ 0.88 was not detected. Our mass maps of 
clusters SPT0517-5430, SPT0509-5342, and SPT0528-5300 have 
peak S/N=3.2, 3.0 and 2.5 respectively. The mass reconstruction 
of SPT0528— 5300, at z ~ 0.7, is shown in the lower panel of 
Fig. Q] There is a second peak (S/N = 3.5) which may warrant 
further investigation; it is 6.5' from the brightest cluster galaxy 
(BCG). It is at 05 h 27 m 46.2 s , -53 h 07 m 57.9 s and we estimate it 
t .obe8 .12±^°i x 10 14 M . On examination of Staniszewski et al. 
(2008, Fig.l), a small decrement in the SZ appears to be present. 
There is no statistically significant difference between the redshift 
distribution of galaxies in a 3' radius around this location compared 
to the redshift distribution of the surrounding 0.4 square degree re- 
gion. 

Figure [2] shows the average convergence in apertures for the 
three detected SPT clusters and for the offset peak. In red (empty 
circles) we show the B-mode which suggests that the shear cata- 
logues are reasonably free of systematics. Figure[3]compares lens- 
ing and optical mass estimates. We see that the SZ-discovered 
detections are amongst the most massive of the clusters ) 3 x 
1O 14 M0. Ml 200 is an estimate of the mass, calibrated b y weak 
lensing measurements for SDSS clusters dReves et alj2008h . within 
a radius r2oo in which the number density is estimated to be 
200/f2 m times the average galaxy number density. This is sub- 
ject to uncertainties in bias, but is claimed to be an unbiased (to 
5%) estimate o f the radius where the mass density is 200 times the 
critical density (Johnston et al. 2 0071). The uncertain ty in Ml 200 is 
estimated to be a factor of 2 ( Menante au et alj2009l) due to the un- 
certainty in extrapolating the scaling relation to higher redshifts and 
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Figure 1. Lensing mass reconstructions for SZ-detected clusters. Contours 
show the projected lensing convergence (mass distribution), at 1%, 1.5%, 
2%, ... 4%. A circle indicates the BCG in each optically selected cluster, an 
X indicates a peak in X-ray emission, and a + denotes the SZ peak. 
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Figure 2. Average convergence within an aperture, h(9) for the three de- 
tected SPT clusters plus the new offset peak found near SPT0528-5300. The 
red (empty circle) shows the B-mode systematic error. 



uncertai n cluster membership. The correlation between the o ptical 
masses dMenanteau et al.l |2009| ; [fvlenanteau & Hughes] l2008h and 
the weak lensing masses gives some justification for using the op- 
tical luminosity as a mass proxy. Interestingly, the most discrepant 
of the SZ clusters has conflicting optical and X-ray mass estimates, 
suggesting that Ml 200 is overestimated. 

We have also calculated the expected Compton y parameter 
and its integral over solid angle, Y = / dfi y(Q) for the clusters. In 
TablesQ]and[2]we show the expected temperature decrement in the 
Rayleigh-Jeans limit (— ATrj) = 2Tcmb{v), averaged within 
r2oo- At 150GHz , the actual decrement is smaller by a factor 0.29 
dCarlstroml20 02, Fig. 2). To compute Y we assume that gas follows 
dar k matter and obtain an a nalytic result. It ignores the gas history 
fcf. iReid & Sperg el (2006)) but our simple model agrees, within 
25%, with th e empirical scaling relations of Mo tl et al] J200 5) and 
Nagai (2006), and is supported by lAtrio-Barandela et alj ( 2008) 
who show from stacked SZ clusters that the baryon profile is con- 
sistent with NFW. Y as a function of M200 is 
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where Da is the angular diameter distance, M200 is the mass 
within an overdensity 820a = 200, a is the total pressure divided by 
the electron pressure. We assume that ions and electrons are in ther- 
mal equilibrium with a Helium mass fraction of 0.24, so a = 1.93. 
We also assume the concentration index 

c s (M 200 ) = 9.59(1 + z)- 1 {M 2m /lQ li h- 1 M & )- - W2 (2) 
dPolag et alj2004h . but note that this is for a as = 0.9 cosmology, 
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Figure 3. Lensing mass A/200 measurement against optical mass Afj J900 
for SZ and optically selected clusters. Filled circles show the SZ-selected 
and optically confirmed clusters, while empty circles denote clusters ob- 
served optically. Lower error bars marked with a triangle signify that the 
lower error reaches zero. X-ray mass estimates are shown as an X on this 
plot for the 3 SZ-detected clusters, joined by an arc to the optical estimate. 
The dashed line is M^ 200 = M20 to guide the eye. The uncertainty in 
Ml is estimated to be a factor 2 iMenanteau et alj20O9l) . 



higher than current estimates jKomatsuetail fcorm Values of Y 
were predicted for clusters with mass detections, with statistical 
errors which are dominated by the error in M^$ a . 

With a uniform prior on masses, we record in Table Q] the 
most likely masses and asymmetric one-sigma errors for the SPT 
SZ clusters. X-ray masses from soft X-ray luminosity using the 
correlations in Reiprich & Bohringer (2002) are shown in TableQ] 
Predicted Y and average y within V2oo are also shown here, along 
with the cluster centre used. To the extent that we can estimate (y) 
from the published maps, our results seem to be about 50% higher 
than observed, but within the errors expected from the mass de- 
terminations. Recent ACT SZ measurements dHincks et af]|2009h 
are in agreement with our predictions for the two low-redshift clus- 
ters which we can compare. Observations of SPT05 17-5430 with 
XMM-Newton yiel d an X-ray mass within rsoo of 6.4 x 10 14 M Q 
dZhang et alj|2006h . Table [2] shows mass measurements and pre- 
dicted Y and average y within r2oo for the optically-detected clus- 
ters. 



5 CONCLUSIONS 

There has been longstanding optimism that SZ selection would be 
among the most favourable ways of detecting clusters for cosmo- 
logical studies, since simulations show that the SZ detection thresh- 
old corresponds very nearly to a threshold in mass. Despite this, 
it was not guaranteed that the first SZ experiments could trace 
mass. In order to demonstrate this we have measured, for the first 
time, weak lensing masses of clusters detected by their Sunyaev- 
Zel'dovich signature. Of the four clusters detec ted by the SPT and 
published recently by St aniszewski et al] d200 8h . we have detected 
three of them, using optical imaging data from the Southern Cos- 
mology Survey. The fourth cluster, at redshift of 0.88, is too distant 
to be detected with these optical data. By fitting NFW profiles we 
have established that the published SZ peaks correspond to mass 
regions, and so we can confirm that the first installment of SZ se- 
lected clusters trace the most massive clusters. 

We have also presented weak lensing mass estimates for other 
clusters detected optically in the SCS. As one might expect the 



published SZ clusters have masses at the upper end of the sample 
10 14 — lO 15 Af0. Furthermore, we have presented analytic predic- 
tions for the integrated Compton Y parameter for the all clusters in 
the sample for future comparison with SZ observations. 
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Table 1. Physical properties of SZ selected clusters in the SCS regions. Redshifts are the photometric redshift of the BCG, with ±1ct limits. The ID is based 
on the position of the BCG. The NFW fit is centred at either the peak in the mass map (map) where available, or the luminosity-weighted centre (lum.). The 
uncertainty in Mr, 200 is estimated to be a factor 2 i Menanteau et al. 2009). Values of 2Tcmb{v) an d Y are predicted from M^jo 1 " usm g Eq. [T] Note that 
ATgz at 150GHz is 0.29 of the penultimate column value. 
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Table 2. As Table 1, but for optically selected clusters in the SCS regions. The NFW fit is centred at the luminosity-weighted centre. 



ID R.A. Dec. 2photo M Lwo M^ " s 2T C M B (y) Y 

Centre of NFW fit (10 14 M Q ) (1O 14 M ) (fiK) ( 10" 5 
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SCSO J232540.2-544430.9 


23:25:32.2 


-54:44:21.5 


0.10±0.02 


2.10 


0.231°;*! 


7+88 


61? 


SCSO J232230.9-541608.3 


23:22:27.2 


-54:16:26.9 


0.12±0.02 


1.62 


o.85i ;H 


97+64 
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SCSO J233000.4-543707.7 


n/a 




0.14±0.02 


1.19 








SCSO J232419.6-552548.9 


23:24:33.6 


-55:26:14.4 


0.18±0.04 


1.19 


< 0.26 


- 


- 


SCSO J233 106.9-555 119.5 


23:31:08.2 


-55:50:56.4 


0.19+.0.04 


0.55 


20+ 69 
u ' zu -0.20 


8lf 


9+19 
—2 


SCSO J233252.9-561454.1 


23:32:51.4 


-56:15:29.8 


0.20±0.03 


1.17 


< 0.09 






SCSO J233215. 5-544211.6 
SCSO J233037. 1-554338.8 
SCSO J233550.6-552820.4 


23:32:17.1 
23:30:34.8 
23:35:46.4 


-54:42:43.1 
-55:43:41.5 
-55:29:21.3 
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23:29:55.8 
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1 1 7+350 
11 ' -110 

92ir 


49+75 

qo+88 
JO -30 
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